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, Abstract. The acceleration of energetic particles in astrophysical shear flows is analyzed. We show 

that in the presence of a non-relativistic gradual velocity shear, power law particle momentum 
distributions f(p) « p-( 3 + a ) may be generated, assuming a momentum-dependent scattering time 
rr ■ T oe p a , with a > 0. We consider possible acceleration sites in astrophysical jets and study the 

conditions for efficient acceleration. It is shown, for example, that in the presence of a gradual shear 
flow and a gyro-dependent particle mean free path, synchrotron radiation losses no longer stop 
the acceleration once it has started to work efficiently. This suggests that shear acceleration may 
naturally account for a second, non-thermal population of energetic particles in addition to a shock- 
accelerated one. The possible relevance of shear acceleration is briefly discussed with reference to 
' the relativistic jet in the quasar 3C 273. 

o ' 

^ INTRODUCTION 

p ' 

Collimated relativistic outflows have been detected in a variety of sources including 
AGNs, microquasars and y-ray bursts (e.g., lfl4l ntX ). Today there is strong evidence 
that many (if not all) of these jets are characterized by a significant shear in their 
velocity field (cf. lEUl for a review). It appears thus very likely that energetic particles 
can be accelerated efficiently by scattering off magnetic inhomogeneities carried within 
such shear flows. Possible implications of such an acceleration process have so far 
been studied by several authors: in the pioneer works on non-relativistic gradual shear 
flows, Berezhko & Krymskii (1981), for example, showed that under suitable conditions 
power law particle momentum distributions may be formed (cf. also |y|]), whereas 
Earl, Jokpii & Morfill (1988) derived the corresponding diffusive particle transport 
equation including shear and inertial effects, assuming the diffusion approximation to 
be valid. The work by Earl et al. was generalized to the relativistic regime by Webb 
(cf. H26L I27L 12811 ) using a mixed-frame approach. Based on those results, Rieger & 
Mannheim (2002) have recently studied the acceleration of particles in rotating and 
shearing flows with application to the relativistic jets in AGN. Particle acceleration in 
non-gradual relativistic shear flows (e.g., in the presence of a relativistic velocity jump), 
on the other hand, has been studied by Ostrowski (| 15, [l6l [l7|]) using Monte Carlo 
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simulations, showing that very flat momentum spectra may be possible. 
In the present contribution we will particularly focus on the conditions required for 
efficient shear acceleration in order to assess whether shear- accelerated particles may 
indeed significantly contribute to the emission at a certain energy band. 




FIGURE 1. Sketch of a two-dimensional, non-relativistic shear velocity profile with the flow directed 
along the z-axis (left). Qualitatively, a particle may be considered as experiencing a gradual shear profile 
if its mean free path A is much smaller than the width Ax of the shear layer and if the flow velocity 
changes smoothly within that layer. In the case where the particle becomes so energetic that its mean 
free path A = cx becomes larger than Ax, it essentially experiences a non-gradual velocity shear, i.e. a 
discontinuous jump in the flow velocity (right). 



BASIC PRINCIPLES OF SHEAR ACCELERATION 

Shear acceleration is based on the idea that energetic particles may gain energy by 
scattering off systematically moving small-scale magnetic field irregularities. These 
irregularities are thought to be embedded in a collisionless shear flow such that their 
velocities correspond to the local flow velocity, i.e. second-order Fermi effects are 
neglected. The scattering process is assumed to occur in such a way that the particles 
are randomized in direction and their energies conserved in the local comoving fluid 
frame. In the presence of a velocity shear, the momentum of a particle travelling across 
the shear changes (with respect to the local fluid frame), so that a net increase may 
occur (e.g., IIOIO . For illustration consider a non-relativistic continuous shear flow with 
velocity profile given by (cf. Fig.lU) 

u = u z (x)e z . (1) 

Let v = (v x ,v y ,v z ) be the velocity vector, m the relativistic mass and p\ the initial 
momentum (relative to local flow frame) of the particle. Within one scattering time 
T (initially assumed to be independent of momentum) the particle travels a distance 
8x = v x % across the shear, so that in the presence of a gradual shear the flow velocity will 
have changed by 8u = 8ue z , with 8u = (dujdx) 8x. Hence the particle's momentum 
relative to the flow becomes P2 = P\ +m8u, i.e. 

p\ = Pi +2m8u pi jZ + m 2 (8u) 2 . (2) 

As the next scattering event preserves the magnitude of the particle momentum (in the 
comoving frame) the particle magnitude will have this value in the local flow frame and 
hence a net increase in momentum may occur with time. Using spherical coordinates 
and averaging over solid angle (assuming an almost isotropic particle distribution), the 
average rate of momentum change and the average rate of momentum dispersion can be 



written as (cf. GjJl) 
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i.e. both depending on the square of the flow velocity gradient. Note that it can be 
shown that a momentum-dependent scattering time obeying a power law of the form 
x oc p a may be accommodated by replacing 4 — > (4 + a) in Eq. ©. Using these results 
we may write down a simple Fokker-Planck transport equation for the phase space 
particle distribution function f(p), assuming a mono-energetic injection of particles 
with momentum pq. Solving for the steady-state with a > one immediately arrives 
at (see JH) 

f(p)~p- {3+a) H(p- P0 ), (5) 

where H(p) is the Heaviside step function. Hence for a mean scattering time scaling with 
the gyro-radius (Bohm case), i.e. x p, Ot = 1, one obtains f(p) p~ 4 , i.e. a power 
law particle number density n(p) p~ 2 which translates into a synchrotron emissivity 



POSSIBLE SHEAR SITES IN ASTROPHYSICAL JETS AND 
ASSOCIATED EFFICIENCIES 

In general we may distinguish at least three possible shear sites in astrophysical jets, cf. 
Rieger & Duffy (2004,2005): 

1. Longitudinal gradual shear: Observationally, there is mounting evidence for an in- 
ternal velocity stratification parallel to the jet axis. In several sources (including 3C353, 
M87 and Mkn 501, cf. O, HE, I22D), for example, at least a two-component velocity 
structure consisting of a fast velocity spine and a slower moving boundary layer is in- 
dicated. Such a velocity structure is also indirectly supported by unification arguments 
for BL Lacs and FR I (|0]) and results from hydrodynamical jet simulations (e.g., Jl|]). 
Moreover, Laing et al. (1999) have argued recently that the intensity and polarization 
systematics in kpc-scale FR I jets are suggestive of a radially (continuously) decreasing 
velocity profile v z (r). In order to estimate the possible particle acceleration efficiency 
in the presence of a longitudinal gradual shear, we may consider a simple realization 
where the flow velocity profile decreases linearly from relativistic to non-relativistic 
speeds over a scale Ar. It can be shown then (see Il2lll ) that the minimum acceleration 
timescale is of the order of 

t 3 (Ar)2 (to 

where %(r) > 1 is the (position-dependent) bulk Lorentz factor of the flow and A the 
particle mean free path. As Ap ro ton A e i ec tron the acceleration of electrons is in general 
much more restricted than the acceleration of protons. In the case where the particle 
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FIGURE 2. Maximum width of velocity shear layer Ar (in parsecs) as a function of the bulk Lorentz 
factor Yb, obtained by requiring the acceleration timescale to be smaller or equal to the synchrotron cooling 
timescale, i.e. f acc < f coo i. The maximum width Ar scales with B~ 3 / 2 and is proportional to the particle's 
rest mass square. The acceleration mechanism is generally much more favourable to protons than to 
electrons. 



mean free path scales with the gyro-radius, i.e. A °= y (Bohm case), the acceleration 
timescale becomes proportional to 1/y, and hence reveals the same dependency on the 
particle Lorentz factor y as the synchrotron (or inverse Compton) cooling timescale. 
Losses are thus no longer able to stop the acceleration process once it has started to work 
efficiently. For illustration, we have plotted in Fig. 2 the maximum width Ar allowed 
for efficient acceleration by balancing ? acc with the synchrotron cooling timescale ? coo i, 
using parameters appropriate for the pc-scale radio jets. For a magnetic field strength 
of B = 0.01 G, for example, a transition width Ar e < 0.003 pc (y, = 10) is required for 
efficient electron acceleration, whereas for protons only Ar p < (m p /m e ) 2 Ar e is needed. 
2. Longitudinal non-gradual shear: When a particle becomes so energetic that its 
mean free path A becomes larger than the width of the velocity transition layer, it 
will essentially experience a non-gradual velocity shear, i.e. a (discontinuous) jump in 
the flow velocity. Ostrowski (1990,1998), for example, has convincingly argued that 
the jet side boundary between the relativistic jet interior and its ambient medium may 
naturally represent a relativistic realization of such a non-gradual velocity shear. Based 
on Monte-Carlo simulations the estimated minimum acceleration timescale is of order 
(cf. Oi). 

v 

^acc~10— provided r„ > Ar, (7) 

c 

where r g denotes the gyro-radius of the particle and Ar is the width of the transition layer. 
Due to the condition r g > Ar an efficient acceleration of electrons appears excluded by 
virtue of their associated rapid radiation losses (e.g., for B = 0.01 Gauss and Ar > 0.001 
pc, i.e. t acc > 10 6 sec, electron Lorentz factors y ~ 2 • 10 10 would be required in order to 



fulfill r g > Ar, implying cooling timescales t coo i ~ 4 ■ 10 2 sec well below the acceleration 
timescale). Efficient acceleration of protons, on the other hand, can be quite possible as 
long as their particle mean free paths remain smaller than the width of the jet. 
3. Transversal gradual shear: Apart from a longitudinal velocity stratification, astro- 
physical jets are also likely to have a significant velocity shear perpendicular to their jet 
axis ('transversal shear'). In particular, several independent arguments suggest that the 
flow velocity field in these jets is characterized by an additional rotational component. 
The strong correlation between the disk luminosity and the bulk kinetic power in the jet 
(e. g.. JI19I1 ). for example, and the observational evidence for a disk-origin of jets (e.g., 



ISLllJ])- suggest that a significant amount of rotational energy of the disk is channeled 
into the jet. Such an internal jet rotation is also implied in theoretical MHD models 
for the origin of jets as magnetized disk winds (e.g., [ 24]). In order to evaluate the ac- 
celeration potential associated with such shear flow, Rieger & Mannheim (2002) have 
recently analyzed the relativistic particle transport in rotating and shearing jets. Based 
on an analytical (mixed-frame) approach for the relativistic Boltzmann equation, using 
a simple (BKG) relaxation scattering term, and assuming validity of the diffusion ap- 
proximation, they considered the acceleration of particle in a cylindrical jet model with 
relativistic outflow velocity v z and different azimuthal rotation profiles. In the case of a 
simple (non-relativistic) Keplerian shear profile they found that local power law spectra 
f(p) 06 p~( 3+a ^ are obtained for t °= p a , a > 0, whereas for more complex rotation 
profiles (e.g., flat rotation) steeper spectra may be possible. 

In order to gain insights into the acceleration efficiency, let us consider a flow veloc- 
ity field with relativistic v z and an azimuthal Keplerian rotation profile of the form 
Q(r) = £2t {r m /r) 3 / 2 , where fij is a constant. The acceleration timescale then scales 
as (see El) 



assuming the flow to be radially confined to r m < r < rj, where r, is the jet radius, 
indicating that efficient particle acceleration generally requires a region with significant 
rotation. Note that for a radially decreasing rotation profile the higher energy emission 
will generally be concentrated closer to the axis (i.e. toward smaller radii). A comparison 
of acceleration and cooling timescales shows that electron acceleration is usually very 
restricted (i.e., only possible for r ~ r in ), whereas proton acceleration appears well 
possible. Note again that for A 7 (Bohm case), the acceleration timescales decreases 
with 7 in the same way as t coo \, suggesting that losses are no longer able to stop the 
acceleration process once it has started to work efficiently. 



Observational and theoretical evidence suggest that astrophysical jets may be character- 
ized by a significant velocity shear. While internal jet rotation, for example, is likely to 
be present at least in the initial parts of the jet, a significant longitudinal velocity shear 
(parallel to the jet axis) prevailing all along the jet might be expected for most powerful 
sources. The acceleration of particles occurring in such shear flows may thus naturally 
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APPLICATIONS 



account for a steady second population of synchrotron-emitting particles, contributing 



to the observed emission in addition to shock- accelerated ones (cf., Ill 111 ). Moreover, 
shear acceleration may also offer a natural explanation for the extended radio and op- 
tical emission observed from several sources. In the case of 3C273, for example, the 
optical spectral index is found to vary only very smoothly along the (large-scale) jet in 
contrast to expectations from simple shock scenarios (cf. I12D). The jet seems to be highly 
relativistic even on kpc scales ( 923 ) so that longitudinal shear acceleration may perhaps 
work efficiently nearly all along the jet. Moreover, there is also evidence for helical bulk 
motion in the large-scale jet ([2]) and internal jet helicity on the VLBI mas-scale and 
below, suggesting that particle acceleration due to internal jet rotation may contribute on 
pc-scale and perhaps even on larger scales. 



Discussions with John Kirk, Alexandre Marcowith and Lukasz Stawarz, and support by a 
Marie-Curie Individual Fellowship (MCIF - 2002 - 00842) are gratefully acknowledged. 
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